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ABSTRACT 

Although many models have been proposed, the physical mechanisms responsible for the formation of low- 
mass brown dwarfs are poorly understood. The multiplicity properties and minimum mass of the brown- 
dwarf mass function provide critical empirical diagnostics of these mechanisms. We present the discovery via 
gravitational microlensing of two very low-mass, very tight binary systems. These binaries have directly and 
precisely measured total system masses of 0.025 Mq and 0.034 Mq, and projected separations of 0.31 AU and 
0.19 AU, making them the lowest-mass and tightest field brown-dwarf binaries known. The discovery of a 
population of such binaries indicates that brown dwarf binaries can robustly form at least down to masses of 
^ 0.02 Mq. Future microlensing surveys will measure a mass-selected sample of brown-dwarf binary systems, 
which can then be directly compared to similar samples of stellar binaries. 
Subject headings: gravitational lensing: micro - binaries: general 



1. INTRODUCTION 

Brown dwarfs (BD) are collapsed objects with masses be- 
low the minimum mass required to fuse hydrogen of ^ 
0.08 Mq. Direct imaging surveys have found that isolated 
BD systems to be f airly ubiquitous in the field as well as in 
young clusters (see iLuhmanI (12012b for a review), with fre- 
quencies rivaling those of their more massive hydrogen-fusing 
stellar brethren. However, it is unclear whether BDs simply 
represent the low-mass extension of the initial collapsed ob- 
ject mass function (IMF), and thus formed via the same pro- 
cesses as stars, or if their formation requires additional physi- 
cal mechanisms. 

The minimum mass of the IMF potentially provides an im- 
portant discriminant between various models of BD forma- 
tion, with very low mass BD binary systems being partic- 
ularly important in this regard. This is because predictions 
for the multiplicity properties of low-mass BDs - frequency, 
mass ratio, and separation as a function of total system mass 
and age - differ significantly depending on the formation sce- 
nario. Nowever, the currently available observational samples 
are strongly influenced by detection biases and selection ef- 



fects. For example, although BD surveys in young stellar as- 
sociations allow for detections of low-mass BD systems, these 
associations are typically fairly distant, making it difficult to 
detect tight BD binaries. Conversely, field BD binaries can be 
resolved to much smaller separations, but low-mass, old field 
BDs are quite faint and thus difficult to detect. As a result, 
the current sample of binary BDs is not only small in number 
but is also substantially incomplete, particularly in the regime 
of low mass and small separation. A further complication is 
that direct mass measurements are available only for a subset 
of tight field BD binaries. Mass estimates of other systems 
must rely on comparison with models, resulting in substantial 
systematic uncertainties. 

Gravitational microlensing is well-suited to fill the gap. 
Microlensing is the astronomical phenomenon wherein the 
brightness of a star is magnified by the bending of light caused 
by the gravity of an intervening object (lens) located between 
the background star (source) and an observer Since this effect 
occurs regardless of the lens brightness , microlensing is suit- 
able to detect faint objects such as BDs (lPaczvnskj|1986h . For 
a lensing event produced by a binary lens with well resolved 
brightness variation of the lensed star, it is possible to pre- 
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cisely measure the physical parameters of the lensing object 
including the mass and distance. Here we report the discovery 
and characterization of two binary BD systems, which both 
have very low mass and tight separation, thus constituting a 
new population. 

2. OBSERVATION 

These BD binaries were discovered in microlens- 
ing events OGLE-2009-BLG-151/MOA-2009-BLG- 
232 and OGLE-201 l-BLG-0420. The events oc- 
curred on stars located in the Galactic Bulge field with 
equatorial and Galactic coordinates (RA.DEC)2oo() = 
(17*'54"22.34S-29°03'20.8"), (/,^)2()oo = (0.88°, -1.70°) 
and (RA,DEC)2()oo = (17^'50™56.18^-29°49'30.2"), 
(l:b)200Q = (359. 84°, -1.45°), respectively. 

OGLE-2009-BLG-151/MOA-2009-BLG-232 was first dis- 
covered by the Optical Gravitational Lensing Experiment 
(OGLE: Udalskj (2003J) group and was independently dis- 
covered b}Mhe_^ficro]^sing_Observation^ Astrophysics 
(MOA: lBOTid7etan (l200lb : ISumi et all (l2003i) ) group in 2009 
observation season. OGLE-201 l-BLG-0420 was detected by 
the OGLE group in the 2011 season. Both events were ad- 
ditionally observed by follow-up observation groups includ- 
ing Microlensing Follow-Up Network (/iFUN: Gould et al] 
(|2006i)). Probing Le nsing Anoma Ues NETwork (PLA NET: 
( iBeauheu et £11120061) ). RoboNet ( Tsapras et all (l2009h ). and 
Microlensing Network for the Detection of Sm all Terrestrial 
Exoplanets (MiNDSTEp: Dominik et al] ( l2010h ). In Table[l] 
we list the survey and follow-up groups along with their tele- 
scope characteristics. Data reductions were carried out using 
photometry codes developed by the individual groups. 

3. ANALYSIS 

Figure 1 displays the light curves of the individual 
events. OGLE-2009-BLG-151/MOA-2009-BLG-232 is char- 
acterized by two strong spikes flanking a "U"-shape trough, 
which is typical for caustic-crossing binary-lens events. Caus- 
tics denote positions on the source plan e where the lensin g 
magnification of a point source diverges (iPetters et al.ll200lT) . 
When a caustic is formed by an astronomical object composed 
of 2 masses, it forms a single or multiple sets of closed curves 
each consisting of concave curves that meet at cusps. When 
a source star crosses the caustic, its brightness is greatly en- 
hanced, causing strong deviation from the smooth and sym- 
metric single-lens light curve. The light curve of OGLE- 
2011-BLG-0420, on the other hand, appears to be smooth and 
symmetric, which are the characteristics of a lensing event 
caused by a single mass. From the fit based on the single 
lens model, however, the light curve exhibits noticeable devi- 
ations near the peak, which indicates the existence of a com- 
panion to the lens. According to the classification scheme of 
binary signatures in lensing light curves set by Ingrosso et al] 
(|2009), the deviations of OGLE-2009-BLG-151/MOA-2009- 
BLG-232 and OGLE-201 l-BLG-0420 are classified as Class 
II and I, respectively. 

With known binary signatures, we conduct binary-lens 
modeling of the observed light curves. For a single lens, the 
light curve is described by 3 parameters: the time of closest 
lens-source approach, to, the lens-source separation (normal- 
ized by the Einstein-ring radius, 9e) at that time, mq, and the 
Einstein time scale, ?£, which represents the time required for 
the source to cross 9e- The Einstein ring denotes the image 
of a source in the event of perfect lens-source alignment, and 
so is used as the length scale of lensing phenomena. Binary 
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Fig. 1 . — Light cui^ves of the binary BD microlensing events OGLE-2009- 
BLG-151/MOA-2009-BLG-232 and OGLE-201 l-BLG-0420. For the MOA 
data of OGLE-2009-BLG-151/MOA-2009-BLG-232, binned data are addi- 
tionally plotted to better show residuals. 
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Fig. 2. — Geometiy of the lens systems. In each panel, the cuspy closed fig- 
ure represents the caustic, the small filled dots are the locations of the binary 
lens components, and the curve with an arrow represents the source trajec- 
tory. Two sets of lens positions and the con'esponding caustics are presented 
at the times marlced in the panel. The circle on the source trajectory repre- 
sents the scale of the source star with respect to the caustic. AH lengths are 
normalized by the Einstein radius con'esponding to the total mass of the lens. 

lenses require three additional parameters: the mass ratio, q, 
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TABLE 1 
Telescopes 



event 


telescope 


OGLE-2009-BLG-151 


MOA, 1.8ra Mt. John, New Zealand 


/MOA-2009-BLG-232 


OGLE, 1.3m Warsaw, Las Campanas, Chile 




^tFUN, 1.3m SMARTS, CeiTO Tololo Inter- American (CTIO), Chile 




^tFUN, 1.0m Mt. Lemmon, USA 




PLANET, l.OmCanopus, Australia 




PLANET, 1.0m South Afiican Astronomical (SAAO), South Africa 




PLANET, 0.6m Perth, Australia 




MiNDSlbp, 1.54m Danish, La Silla, Chile 




RoboNet, 2.0m Faulkes North Telescope (FTN), Hawaii, USA 




RoboNet, 2.0m Faulkes South Telescope (FTS), Hawaii, USA 




RoboNet, 2.0m Liverpool Telescope (LT), Canary Islands, Spain 


OGLE-2011-BLG-0420 


OGLE, 1.3m Warsaw, Las Campanas, Chile 




MOA, 0.6m B&C, Mt. John, New Zealand 




^iFUN, 1.3m SMARTS, CTIO, Chile 




/iFUN, 1.0m Wise, Israel 




PLANET, 1.0m SAAO, South Africa 




PLANET, l.OmCanopus, Australia 




MiNDSTEp, 1.54m Danish, La Silla, Chile 




RoboNet, 2.0m LT, Canary Islands, Spain 



the projected separation (normalized by 6^) between the bi- 
nary components, s, and the angle between the source trajec- 
tory and the binary axis, a (source trajectory angle). 

In addition to the basic lensing parameters, it is often 
needed to include additional parameters to precisely describe 
subtle light-curve features caused by various second-order 
effects. For both events, the lensing-induced magnification 
lasted for several months, which comprises a significant frac- 
tion of the Earth's orbital period around the Sun (1 year). 
Then, the apparent motion of the source with respect to the 
lens deviates from rectilinear (Gould 1992) due to the change 
of the observer's position caused by the Earth's orbital mo- 
tion. This parallax effect causes long-term deviation in lens- 
ing light curves. Consideration of the parallax effect requires 
2 additional parameters of tte.w and tte.e, which are the two 
components of the lens parallax vector tte, projected on the 
sky along the north and east equatorial coordinates, respec- 
tively. The orbital motion of the lens also affects lensing 
light curves. The lens orbital motion causes the projected bi- 
nary separation and the source trajectory angle to change over 
the course of a lensing event. These require two additional 
lensing parameters of the change rates of the binary separa- 
tion, ds/dt, and the source trajectory angle, da/dt. Finally, 
finite-source effects become important whenever the magnifi- 
cation varies very rapidly with the change of the source po- 
sition, so that different parts of the source are magnified by 
different amounts. Such a rapid magnification variation oc- 
curs near caustics and thus finite-source effects are important 
for binary-lens events involved with caustic crossings or ap- 
proaches. This requires one more parameter, the normalized 
source radius = 6'h./6'e, where 9^. is the angular source ra- 
dius. Measuring the deviation caused by the parallax and 
finite-source effects is important to determine the physical pa- 
rameters of the lens. By measuring the finite-source effect, the 
Einstein radius is determined by 6-e = 0*1 once the source 
radius is known. With the measured lens parallax and the 
Einstein radius, the mass and distance to the lens are deter- 
mined as Mtot = fe/(K7rE) and Z?l = A U /(tte^e + tts ), respec- 
tively (Gould" 1992( lGould et al.ll2006h . Here k = 4G/(c2aU), 
AU is an Astronomical Unit, tts = AU/Z)s, and Ds '--^ 8 kilo- 
parsec is the source distance. 

We model the observed light curves by minimizing x~ in the 



TABLE 2 
Best fit lensing parameters 



parameters 


OGLE-2009-BLG-151 
/MOA-2009-BLG-232 


OGLE-2011-BLG-0420 


X-/dof 


3040.9/3032 


5410.3/5439 


?o (HID-) 


4999.680 ± 0.061 


5766.110 ±0.001 


"0 


-0.217 ± 0.004 


-0.030 ± 0.001 


;e (days) 


27.95 ±0.11 


35.22 ± 0.08 




1.135 ±0.004 


0.289 ± 0.002 


1 


0.419 ± 0.006 


0.377 ± 0.009 


a 


-1.049 ± 0.004 


-2.383 ± 0.002 


p* (10-2) 


1.06 ±0.01 


4.88 ± 0.01 




-3.33 ±0.11 


-1.15 ±0.05 




-0.91 ±0.11 


0.19 ±0.01 


ds/dt (yr"') 


1.55 ±0.13 


-2.59 ± 0.07 


da/dt (yr"') 


0.69 ± 0.06 


6.88 ± 0.20 


slope (mag yr"' ) 


0.0052±0.0008 





Note. — HJD'=HJD-2450000. 



parameter space. We investigate the existence of possible de- 
generate solutions because it is known that light curves result- 
ing from different combi nations of lensing parameters often 
result in a simil ar shape tGriest & Safazadehlll998l DominikI 
IT999HAnll2005h . In modeling finite -source effects, we addi- 
tionally consider the limb-darkening variation of the source 
star surface (Witt 1995) by modeling the surface profile as a 
standard linear law. For minimization, we use the Markov 
Chain Monte Carlo method. Photometric errors of the indi- 
vidual data sets are rescaled so that per degree of freedom 
becomes unity for each data set. We eliminate data points 
with large errors and those lying beyond 3 a from the best-fit 
model to minimize their effect on modeling. 

Table|2]gives the solutions of the lensing parameters found 
from modeling. In Figure 2, we also present the geometry 
of the lens system where the source trajectory with respect 
to the positions of the binary lens components and the caus- 
tic are shown. For OGLE-2009-BLG-151/MOA-2009-BLG- 
232, we find that the two strong spikes were produced by 
the source crossings of a big caustic formed by a binary lens 
with the projected separation between the lens components 
(s ^ 1.14) being similar to the Einstein radius of the lens. We 
find that including the second-order effects of lens parallax 
and orbital motions improves the fit by Ax^ =213. OGLE- 
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TABLE 3 
Source Star Properties 



quantity 


OGLE-2009-BLG-151 
/MOA-2009-BLG-232 


OGLE-2011-BLG-0420 


(V-/)o 


1.352 


1.611 


k 


14.490 


13.091 


Ot (pas) 


7.57 ± 0.66 


15.94 ± 1.38 


stellar type 


K giant 


K giant 


Be (mas) 


0.71 ±0.01 


0.33 ± 0.03 



201 l-BLG-0420 is also a caustic-crossing event, but the pro- 
jected binary separation (s ^ 0.29) is substantially smaller 
than the Einstein radius. For such a close binary lens, the 
caustic is small. For OGLE-2011-BLG-0420, the caustic is 
so small that the source size is similar to that of the caus- 
tic. Hence, the lensing magnification is greatly attenuated 
by the severe finite-source effect and thus the deviation dur- 
ing the caustic crossings is weak. We find that there exists 
an alternative solution with > 1 caused by the well-known 
close/wide binary degeneracy, but the degeneracy is resolved 
with Ax^ = 27. The parallax and lens orbital effects are 
also clearly measured with Ax^ = 403. From 10 years of 
OGLE data, we find that OGLE-2011-BLG-0420S (source 
star) is extremely stable, but OGLE-2009-BLG-151/MOA- 
2009-BLG-232S exhibits irregular < 1 % variations, typically 
on time scales of a few hundred days. Because these can affect 
the parallax measurement, we restrict the modeling to fo ± 300 
days to minimize the impact of variations while still retain- 
ing enough baseline to ensure a stable fit. We also include 
a "slope" parameter for the source flux to account for the re- 
maining variability. We find only slight differences in final re- 
sults if we repeat this procedure with longer baselines. There- 
fore, it is unlikely but not impossible that source variability 
affects the OGLE-2009-BLG-151/MOA-2009-BLG-232 par- 
allax measurement. By contrast, the results for OGLE-2011- 
BLG-0420 are very robust. 

4. PHYSICAL PARAMETERS 
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Fig. 3. — Locations of the lensed stars of the individual events on the color- 
magnitude diagrams. 

Since and tte are well measured (Table |2]l, it is pos- 
sible to determine Mtot and Dl for both systems. The only 
missing ingredient is the angular source radius 9^,, needed to 
find 9e = This is determined from the de-reddened 

color and brightness of the source star For the calibration 
of the color and brightness, we use the cent roid of bulge gi- 
ant clump as a reference (lYoo et alj 120041) because its de- 
reddened brightness /q.c = 14.45 at the Galactocentric distance 
dNataf et al..201Z) and the color (y-/)o,c = 1 -06 (.Bensbv et ah 



TABLE 4 
Physical quantities 



quantity 


OGLE-2009-BLG-151 
/MOA-2009-BLG-232 


OGLE-2011-BLG-0420 


Mtot (Mq) 


0.025 ± 0.001 


0.034 ± 0.002 


Ml (Mq) 


0.018 ±0.001 


0.025 ± 0.001 


M2 (Mq) 


0.0075 ± 0.0003 


0.0094 ± 0.0005 


Dl (kpc) 


0.39 ± 0.01 


1.99 ±0.08 


Jj_ (AU) 


0.31 ±0.01 


0.19 ±0.01 



l20Tlh are known. We then translate V -I into V -K color 
by using the relation (Bessell & Brett 1988) and then find 9^, 
using the relation be tween the V -K and the angular radius 
(iKervella etalj|200l . In Tabled we list the measured de- 
reddened colors (V-I)o, magnitudes /q, angular radii, types 
of the source stars and the measured Einstein radii of the indi- 
vidual events. Figure|3]shows the locations of the lensed stars 
of the individual events on the color-magnitude diagrams. 

The derived physical quantities for the OGLE-2009-BLG- 
151/MOA-2009-BLG-232L and OGLE-2011-BLG-0420L 
binaries are listed in Table |4] Here the letter "L" at the 
end of each event indicates the lens of the event. The total 
system masses are Mtot = (0.025 ± 0.001) Mq and (0.034 ± 
0.002) Mq, respectively, well below the hydrogen-burning 
limit. The projected separations and mass ratios are d± = 
(0.31 ±0.01) AU and q = 0.419 ±0.006 for OGLE-2009- 
BLG-151/MOA-2009-BLG-232L, and dj_ = (0.19 ±0.01) AU 
and q = 0.377 ± 0.009 for OGLE-201 1-BLG-0420L. It is 
worth emphasizing the high precisions (< 10%) with which 
the total system masses and individual component masses are 
determined. 

Figure m compares OGLE-2009-BLG-151/MOA-2009- 
BLG-232L and OGLE-201 1-BLG-0420L to a sample of 
low-mass binaries in the field and in young a ssociations from 
iFahertv et a P (201 1). Basri & Martini (ll999').'Burgasser et al J 
(120081) iBurgas ser e t al. (2012|), and Lane et al. (2001). 
The only known BD binaries with comparable total 
masses are Oph 16225-240515 with Mtot ^ 0.032 Mq 
dlavawardhana & Ivanov 2006) and 2MASSJ1207334- 
393254 ( Chauvin etal.i ,2004) with Mtot 0.028 Mq. 
However, these two systems are both young (5 Myr and 
8 Myr, respectively), and have very wide separations of 
hundreds of AU. Indeed, OGLE-2009-BLG-151/MOA-2009- 
BLG-232L and OGLE-201 1-BLG-0420L are the tightest 
known BD binaries with substantially lower mass than 
previously known field BD binaries. Both systems have mass 
ratios of ^ 0.4, apparently consistent with the trend found 
from old field BDs, which tend t o have a prefer e nce f or 
larger mass ratios (see Figure |3] of IBurgasser et al.l (l2007t ). 
although it is important to stress that the selection effects in 
m icrolensing and direct imaging surveys are very different. 

'Burgasse r et al.l (120071) suggested that field low-mass bi- 
naries with Mtot = 0.05 - 0.2 Mq may exhibit an empirical 
lower limit to their binding energies of Gm i ma/a 2 .5 x 1 0"*^ 
erg (see Figure E). Although OGLE-2009-BLG-151/MOA- 
2009-BLG-232L and OGLE-201 1-BLG-0420L are substan- 
tially lower in mass than these BD binaries, they are also 
considerably tighter Therefore, with binding energies of 
'--^ 7 X 10"*^ erg and 2 x lO'*-' erg, they are consistent with the 
extrapolation of the minimum binding energy limit down to 
total system masses of Mtot ^ 0.02 Mq. 

Although we are unable to provide an estimate of the space 
density of such tight, low-mass brown dwarf binaries, nor 
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Fig. 4. — Projected separation versus total system mass for a compilation 
of binaries. Grey circles indicate old field binaries, whereas blue squares 
indicate young (< 500 Myr) systems. The size of the symbols is proportional 
to the square root of the mass ratio. The red stars are the two tight, low-mass 
binary BDs discussed here, which have mass ratios of ~ 0.4. The dashed line 
shows a binding energy of 2 X 10*^ erg, assuming a mass ratio of 1. 
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Fig. 5. — Binding energy (Gmim2 /a) versus total system mass for the same 
binaries as shown in Fig. |4] Symbols and Une are the same. 

even an estimate of their frequency relative to more massive 
stellar binaries, the discovery of two systems among the rela- 
tively small sample of binary lensing events with precise mass 
estimates strongly suggests that very low-mass, very tight BD 
binaries are not rare. Thus these detections herald a much 
larger population of such systems. We can therefore conclude 
that BD binaries can robustly form at least down to system 
masses of ^ 0.02 Mq, providing a strong constraint for for- 
mation models. 



5. DISCUSSION 

The discoveries of the binary BDs reported in this paper 
demonstrate the importance of microlensing in BD studies. 
The microlensing method has various advantages. First, it en- 
ables to detect faint old populations of BDs that could not 
be studied by the conventional method of imaging surveys 
and the sensitivit y extends down to planetary mass objects 
(ISumi et al.l201 V ). It also allows one to detect BDs distributed 
throughout the Galaxy. Therefore, microlensing enables to 
study BDs based on a sample that is not biased by the bright- 
ness and distance. Second, in many cases of microlensing 
BDs, it is possible to precisely measure the mass, which is 
not only the most fundamental physical parameter but also 
a quantity enabling to unambiguously distinguish BDs from 
other low mass populations such as low mass stars. While 
mass measurements by the conventional method require long- 
term and multiple stage observation of imaging, astrometry, 
and spectroscopy by using space-borne or very large ground- 
based telescopes, microlensing requires simple photometry by 
using 1 m class telescopes. Despite the observational sim- 
plicity, the mass can be measured with uncertainties equiva- 
lent to or smaller than those of the measurement by conven- 
tional methods. Finally, microlensing can expand the ranges 
of masses and separations in the binary BD sample that is 
incomplete below ^0.1 Mq in mass and ^ 3 AU in separa- 
tion. Microlensing sensitivity to binary objects peaks when 
the separation is of order the Einstein radius. Considering 
that the Einstein radius corresponding to a typical binary BD 
is < 1 AU, microlensing method will make it possible to study 
binary BDs with small separations. 

The number of microlensing BDs is expected to increase in 
the future with the upgrade of instruments in the existing sur- 
vey experiments and the advent of new surveys. The OGLE 
group recently upgraded its camera with a wider field of view 
to significantly increase the observational cadence. The Korea 
Microlensing Telescope Network (KMTNet), now being con- 
structed, will achieve 10 minute sampling of all lensing events 
by using a network of 1 .6 m telescopes on three different con- 
tinents in the Southern hemisphere with wide-field cameras. 
Furthermore, there are planned lensing surveys in space in- 
cluding EUCLID and WFIRST. With the increase of the mi- 
crolensing event rate combined with the improved precision 
of observation, microlensing will become a major method to 
study BDs. 
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